 Matlab 1 Lansarea aplicatiei Matlab 5 1 1 Fereastra aplicatiei Matlab 5 1 2 Fereastra editorului Matlab 8 2 Elementele limbajului Matlab 12 2 1 Constante 12 2 2 Variabile 12 2 3 Operatori aritmetici Expresii aritmetice 12 2 4 Operatori relationali Expresii relationale 14 2 5 Operatori booleeni Expresii booleene 14 2 5 1 Operatori booleeni ce conecteaza expresii relationale 14 2 5 2 Operatori booleeni ce actioneaza asupra scalarilor si elementelor matricelor 15 2 6 Prioritatile si asociativitatea operatorilor 15 2 7 Operatii cu siruri de caractere 16 2 8 Vectori si matrice 17 2 8 1 Operatii cu vectori si matrice 17 2 8 2 Generarea vectorilor 18 2 8 3 Lucrul cu matrice împartite în blocuri 19 2 8 4 Concatenarea matricelor 19 2 8 5 Functii ce creaza matrice speciale 20 2 8 6 Modificarea dimensiunilor matricelor 21 2 9 Structuri 21 2 10 Functiile load si save 22 2 11 Introducerea datelor de la tastatura 22 2 12 Functii pentru lucrul cu siruri de caractere 22 2 13 Formatarea rezultatelor 24 2 14 Celule Tablouri de celule 24 3 Instructiuni Matlab 28 3 1 Instructiunea for 28 3 2 Instructiunea while 29 3 3 Instructiunile continue si break 30 3 4 Instructiunea if 30 3 5 Instructiunea switch 31 4 Functii 32 4 1 Definirea functiilor 32 4 2 Instructiunea global 32 4 3 Variabile handle de functii 33 4 4 Vectorizarea functiilor 37 4 5 Aplicarea functiilor asupra tablourilor de celule 37 5 Reprezentari grafice 39 5 1 Reprezentarea curbelor plane 39 5 2 Editarea unei figuri 46 5 3 Modificarea interactiva a proprietatilor obiectelor grafice 47 5 3 1 Modificarea proprietatilor axelor 48 5 3 2 Modificarea proprietatilor curbelor 50 5 4 Modificarea programata a proprietatilor obiectelor grafice 54 5 4 1 Obiectul Root 55 5 4 2 Obiectul Figure 55 5 4 3 Obiectul Axes 56 5 4 4 Obiectul Line 57 5 4 5 Obiectul Text 57 5 4 6 Modificarea proprietatilor obiectelor grafice 61 6 Transformata Fourier discreta 66 7 Integrarea ecuatiilor diferentiale ordinare 69 8 Filtre numerice 73 9 Rezolvarea ecuatiilor algebrice neliniare 75 10 Interpolarea functiilor 78 11 Operatii cu polinoame 79 12 Interfete grafice 80 12 1 Construirea unei interfete grafice cu aplicatia GUIDE 80 12 1 1 Descrierea aplicatiei GUIDE 80 12 1 2 Crearea interfetei grafice 85 12 1 3 Alinierea componentelor 87 12 1 4 Modificarea proprietatilor obiectelor grafice 88 12 2 Programarea interfetei grafice 94 12 2 1 Organizarea interfetei grafice 94 12 2 2 Programarea interfetei grafice 100 12 2 3 Adaugarea unor grafice la interfata 102 12 2 4 Modificarea proprietatilor obiectelor grafice la rularea programului 107 12 2 5 Pozitionarea unui punct pe o curba 112 12 2 6 Introducerea datelor cu casete de dialog 119 13 Calcul simbolic 123 13 1 Definirea variabilelor si functiilor simbolice 123 13 2 Substitutii 125 13 3 Reprezentarea simbolica a numerelor 127 13 4 Definirea de variabile simbolice reale si complexe 128 13 5 Matrice simbolice 129 13 6 Crearea de functii abstracte 132 13 7 Reprezentarea grafica a functiilor 132 14 Functii pentru calcule simbolice 136 14 1 Derivarea 136 14 2 Limite de functii 138 14 2 1 Calculul limitelor de functii 138 14 2 2 Calculul limitelor laterale 139 14 3 Integrarea 139 14 3 1 Calculul integralei nedefinite 139 14 3 2 Calculul integralei definite 140 14 3 3 Calculul integralelor ce depind de un parametru 141 14 4 Calculul sumelor 142 14 5 Dezvoltarea în serie Taylor 143 14 6 Simplificari si substitutii 144 15 Operatii cu polinoame simbolice 146 16 Algebra liniara simbolica 149 17 Rezolvarea simbolica a ecuatiilor algebrice 151 18 Rezolvarea simbolica a ecuatiilor diferentiale ordinare 154 19 Transformari integrale 156 19 1 Transformata Fourier 156 19 2 Transformata Laplace 157 19 3 Transformata Z 158 1 Lansarea aplicatiei Matlab Lansarea aplicatiei Matlab se face: * cu un clic pe icoana aplicatiei sau, * cu un clic pe butonul Start si apoi pe optiunea Programs; în lista de aplicatii afisate se executa un clic pe aplicatia Matlab 1 1 Fereastra aplicatiei Matlab Fereastra aplicatiei Matlab este cea din Figura 1 Ea cuprinde, de regula, doua ferestre: * fereastra de comenzi (Command Window) unde se introduc instructiunile Matlab, * fereastra cu instructiunile introduse anterior, (Command History) Figura 1 Fereastra aplicatiei Matlab Utilitatea fereastrei cu instructiunile introduse anterior, (Command History) este aceea ca, la un clic dublu pe o instructiune din aceasta fereastra, ea este copiata în fereastra de comenzi (Command Window), si executata Este posibil de a afisa si alte ferestre: * directorul de lucru Matlab, (Current Directory), ce afisaza continutul directorului current al aplicatiei, * fereastra Help, * fereastra Workspace, ce afisaza variabilele din program Meniurile aplicatiei Matlab au, în general, optiunile cunoscute Meniul File, aratat în Figura 2, are optiunile: * New, crearea unui fisier nou, * Open, pentru a deschide un fisier existent * Close Command History, închide fereastra Command History, * Page Setup, alege optiunile de tiparire a ferestrei active, * Print, tipareste ferestra activa, * Exit Matlab, închide aplicatia Matlab Sunt afisate de asemenea ultimele patru fisiere deschise Aplicatia Matlab poate fi închisa si cu comanda quit sau exit în fereastra de comenzi Meniul File Meniul Edit Figura 2 Meniurile File si Edit ale aplicatiei Matlab Meniul Edit contine optiunile cunoscute: * Undo, anuleaza ultima comanda, * Redo, executa ultima comanda, * Cut, muta textul selectat în Clipboard, * Copy, copiaza textul selectat în Clipboard, * Paste, copiaza textul din Clipboard, * Select All, selecteaza tot textul ferestrei active, * Find, permite cautarea în text a unui sir de caractere, * Clear Command Window, sterge continutul ferestrei de comenzi, * Clear Command History, sterge continutul ferestrei cu comenzile introduse Meniul Vindow, afisat în Figura 3, are optiuni de afisare sau ascundere a diverselor ferestre: * Command Window, afisaza / ascunde ferestra de comenzi, * Command History, afisaza / ascunde ferestra cu comenzile precedente, * Current Directory, afisaza / ascunde directorul de lucru, * Workspace, afisaza / ascunde spatiile de lucru, Figura 3 Meniul Window Fiecare fereastra are un buton de închidere , un buton undock, , care muta fereastra respectiva în afara Ferestrei Matlab si un buton dock, , muta fereastra înapoi în cadrul ferestrei Matlab Meniul Desktop din Figura 4 permite organizarea ferestrei aplicatiei Figura 4 Meniul Desktop Putem afisa doar ferestrele dorite, cu optiunile Command Window, , Help Submeniul Desktop Layout permite de asemenea alegerea ferestrelor dorite si aranjarea lor Bara de instrumente a aplicatiei Matlab este aratata în Figura 5 O parte din butoane corespund optiunilor meniurilor File si Edit Figura 5 Bara de instrumente a aplicatiei Matlab Semnificatia butoanelor este cea din Tabelul 1 Buton Semnificatie Crearea unui nou fisier m Deschide un fisier Cut Copy Paste Undo Redo Lansare Simulink Lansare GUIDE Help Selectarea directorului curent Alegerea unui nou director curent Tabelul 1 Butoanele barei de instrumente a aplicatiei Matlab 1 2 Fereastra editorului Matlab Functiile limbajului Matlab sunt scrise în fisiere cu extensia m cu ajutorul editorului Matlab Ferestra editorului este cea din Figura 6 Figura 6 Fereastra editorului Matlab Meniul File al editorului are optiunile cunoscute, de a deschide un fisier nou sau unul existent, de a salva si închide fisiere, etc Sunt afisate si ultimele patru fisiere deschise Meniul File este aratat în Tabelul 2a Meniul Edit are optiunile Cut, Copy, Paste, Undo, Redo, Find, Replace, etc Meniul Edit este aratat în Tabelul 2b a) meniul File b) meniul Edit Tabelul 2 Meniurile File si Edit ale editorului Matlab Meniul Text din Figura 7a are optiuni pentru comentarea si decomentarea textului, indentarea programului, etc Un grup de instructiuni se poate comenta sau decomenta prin selectarea si clic pe optiunile Comment respectiv Uncomment Meniul Debug din Figura 7b permite executarea instructiunilor din fisier cu optiunea Save File and Run a) meniul Text b) meniul Debug Figura 7 Meniurile Text si Debug ale editorului Matlab Comenzi utile pentru lucrul cu editorul sunt: * clc strege ecranul, * format formateaza rezultatele, vezi paragraful Formatarea rezultatelor, * dir listeaza fisierele din directorul curent, * pwd afisaza numele directorului curent 2 Elementele limbajului Matlab 2 1 Constante Constantele utilizate în Matlab sunt numerele reale, numerele complexe si sirurile de caractere Numerele reale se scriu cu partea subunitara separata de cea întreaga de punct Exemple 7 5, -4 0, -4 Numerele reale se pot scrie cu exponent, de forma e±n sau d±n Exemple 1 23e-2, -2 5e+2, -2 5e2 Numere complexe Sunt predefinite constantele complexe i si j cu Numarul complex a+bi se scrie a+bi, a+bj Se poate defini o variabila complexa,de exemplu , cu instructiunea s=sqrt(-1) In acest caz numarul complex a+bi se scrie a+b*s In Matlab exista câteva constante predefinite: pi, eps, NaN, Inf, realmin, realmax Semnificatia lor este cea din Tabelul 3 Constanta Semnificatia pi eps precizia calculelor, 2-52 Inf NaN (not-a-number) este rezultatul operatiilor 0/0, -, / realmin cel mai mic numar real reprezentabil, 2-1022 realmax cel mai mare numar real reprezentabil, 21023 i, j unitatea imaginara Tabelul 3 Constante predefinite în Matlab 2 2 Variabile Numele de variabile sunt formate din litere si cifre si încep cu o litera Variabilele sunt create automat la întâlnirea unui nou nume Variabilele sunt initializate cu instructiuni de forma ident = expresie Mentionam ca literele mici sunt diferite de literele mari De exemplu, variabila x este diferita de variabila X Pentru a afisa valoarea unei variabila se tasteaza numele ei 2 3 Operatori aritmetici Expresii aritmetice Operatorii aritmetici utilizati în Matlab sunt: +, -, *, /, \, ^ Operatori +, -, *, /, sunt cei patru operatori aritmetici cunoscuti Operatia a \ b se defineste ca Operatorul ^ este ridicarea la putere Expresia a^ b este Functiile matematice uzuale sunt cele de mai jos: sin asin sqrt exp abs floor cos acos log ceil tan atan log10 round atan2 Functia abs() calculeaza modulul unui numar real sau complex Functia atan2(x, y) calculeaza valoarea atan(x / y); Prioritatile operatorilor sunt cele cunoscute Operatorii + si - unari au cea mai mare prioritate, urmati de operatorul ^, urmat de urmati de operatorii *, /, \, iar operatorii + si - binari au cea mai mica prioritate Operatorii unari sunt asociativi la dreapta, operatorii binari sunt asociativi este la stânga De exemplu, expresia 2^2^3 = (2^2)^3 si are rezultatul 64 Functia round() rotunjeste la cel mai apropiat întreg, functia floor() rotunjeste in jos, functia ceil() rotunjeste în sus la cel mai apropiat întreg Exemple de expresii aritmetice si valorile lor sunt aratate în Tabelul 4 Expresie aritmetica Valoare 2 ^ 2 5 5 6569 4 / 3 1 3333 4 \ 3 0 75 8 / 4 / 2 = (8 / 4) / 2 1 3 3 ^ 2 ^ 1 5=(3 3 ^ 2) ^ 1 5 35 9370 -2^-3 -0 125 floor(1 9) 1 floor(-1 9) -2 ceil(1 1) 2 ceil(-1 1) -1 round(1 1) 1 round(1 9) 2 Tabelul 4 Expresii aritmetice Pentru lucrul cu numere complexe sunt utile functiile angle real imag conj complex Functiile real() si imag() calculeaza partea reala si imaginara ale unui numar complex, conj() calculeaza conjugatul unui numar complex Functia angle() calculeaza faza în radiani a unui numar complex Modului numarului complex se calculeaza cu functia abs() Faza se poate calcula si cu functiiile atan() si atan2() Pentru numarul complex z = x + yi, faza se calculeaza ca atan(imag(z) / real(z)) sau atan2(imag(z), real(z)) Functia complex(x, y) creaza numarul complex x + iy din numerele reale x si y Fie numarul complex z = 2+3i Operatii cu acest numar complex sunt aratate în Tabelul 5 Expresia aritmetica Valoare real(z) 2 imag(z) 3 abs(z) 3 6056 angle(z) 0 9828 atan(imag(z)/real(z)) 0 9828 atan2(imag(z), real(z)) 0 9828 Tabelul 5 Operatii cu numere complexe In cadrul expresiilor aritmetice se pot utiliza paranteze rotunde, ( si ), pentru a modifica prioritatile operatorilor Pentru a suprima afisarea rezultatelor, la sfârsitul unei instructiuni se scrie operatorul; Continuarea unei instructiuni pe mai multe linii se face cu operatorul 2 4 Operatori relationali Expresii relationale Operatorii relationali din Matlab sunt: , >=, = =, ~= sau lt, le, gt, ge, eq, ne O expresie relationala este formata din doua expresii aritmetice separate de un operator relational Rezultatul unei expresii relationale este true (adevarat) sau false (fals) ce au valorile 1 si repectiv 0 Prioritatea oricarui operator aritmetic este mai mare decât cea a operatorilor relationali Exemple de expresii relationale si valorile lor sunt prezentate în Tabelul 6 Expresie relationala Valoare 2 ~= 2 0 2 5 && b 0 || b + 1 > 2 1 ~(1 > -1) 0 Tabelul 7 Expresii booleene In scrierea expresiilor booleene sunt utile teoremele lui De Morgan: ~(a || b) = ~a && ~b ~(a && b) = ~a || ~b 2 5 2 Operatori booleeni ce actioneaza asupra scalarilor si elementelor matricelor Acesti operatori sunt &, |, si ~ Ei reprezinta operatorii and (si), or (sau), si respectiv not (nu) Orice numar diferit de 0 este considerat ca reprezinta valoarea adevarat, orice numar 0 reprezinta valoarea fals Operatorul ~ este unar si are aceeasi prioritate ca si operatorii + si - unari Operatorii & si | sunt operatori binari Prioritatile lor sunt respectiv & si |, si mai mici decât cele ale operatorilor aritmetici Exemple de expresii booleene sunt prezentate în Tabelul 8 Expresie booleana Valoare ~1 0 ~5 2 0 ~-7 4 0 ~0 1 a & ~ a 0 a | ~a 1 ~1 > -1 1 ~2 > 1 0 Tabelul 8 Expresii booleene In ultimele exemple se evalueaza mai întâi expresia ~1, deoarece operatorul ~ are prioritate mai mare ca >, vezi tabelul cu prioritatile operatorilor de mai jos Exista functia xor cu doua argumente, ce corespunde operatorului xor (sau exclusiv) Functia xor are valoarea true (adevarat) daca un singur argument are valoarea true (adevarat) De exemplu, expresia xor(3, 4) are valoarea zero 2 6 Prioritatile si asociativitatea operatorilor Prioritatile si asociativitatea operatorilor din Matlab sunt cele aratate în Tabelul 9 Operatorii sunt scrisi în ordinea prioritatii Operator Asociativitate +unar, - unar, ~ dreapta ^ stânga *, /, \ stânga +binar, - binar stânga , >=, = =, ~= stânga & stânga | stânga && stânga || stânga Tabelul 9 Prioritatile si asociativitatea operatorilor 2 7 Operatii cu siruri de caractere Sirurile de caractere se scriu între ' (apostrofuri) De exemplu, sirul abc se scrie 'abc' Daca în sir apare caracterul ' acesta se scrie ca '' (doua apostrofuri) De exemplu, sirul a'b se scrie 'a''b' Un sir de caractere este un vector ce contine codurile numerice ale caracterelor Caracterele se reprezinta în Matlab în codul ASCII Indicele primului element din sir este 1 Numarul de caractere se obtine cu functia length() sau cu functia size() De exemplu length('abc') are rezultatul 3, iar size('abc') are ca rezultat matricea Operatorul de concatenare a sirurilor este [] De exemplu ['abcd' '+-*/' '123'] are ca rezultat sirul abcd+-*/123 Câteva functiile pentru lucrul cu caractere sunt urmatoarele: * functia strcat() cu forma strcat(s1, s2, , sn) concateneaza sirurile s1, s2, , sn, * functia strcmp() cu forma strcmp(s1, s2) compara doua siruri si are rezultatul 1 (true) daca sirurile sunt egale si 0 (false) în caz contrar Sirurile sunt egale când contine aceleasi caractere si au aceeasi lungime, * functiile isletter() si isspace() cu formele isletter(s) isspace(s) au ca argument sirul s si ca rezultat un vector de aceeasi dimensiune cu sirul, ce contine 1 pentru caracterele litere / spatii si 0 în caz contrar Caracterele spatii sunt '\t', '\n' si ' ', * functia ischar() cu forma ischar(s) are ca rezultat 1 daca argumentul s este un sir si 0 în caz contrar, * functia findstr() cu forma findstr(s1, s2) cauta sirul cel mai scurt în cel mai lung Rezultatul functiei este un vector cu indicii fiecarei aparitii a unui sir în celalalt * functiile lower() si upper() cu formele lower(s) upper(s) au ca argument sirul s si ca rezultat sirul convertit în litere mici / mari Exemplu Functia findstr('abcdab', 'ab') are ca rezultat vectorul cu indicii aparitiei sirului 'ab' în sirul 'abcdab' 1 5 Exemplu Fie sirurile str1 = 'sir1' ; str2 = 'sir2' Expresia strcmp(str1, str2) are valoarea zero deoarece sirurile nu sunt egale Exemplu Expresia strcat ('xyz', 'a+x') are ca rezultat sirul xyza+x 2 8 Vectori si matrice 2 8 1 Operatii cu vectori si matrice Definirea unui vector sau a unei matrice se face scriind elementele între paranteze drepte, [ ] Operatorul [ ] se va numi operator de concatenare Exemplu de vector [-1 3, sqrt(3), (1+2)*4/5] Caracterele , si spatiu separa coloanele Caracterele ; si CR (cariage return) separa liniile De exemplu, A = defineste o matrice cu 3 linii si 3 coloane Prima linie este separata de cea de a doua prin ;, a doua linie este separata de a treia prin CR Indicii vectorilor si matricelor încep de la 1 Elementul al matricei a se scrie a(i, j) Operatiile cu matrice se efectueaza cu operatorii +, -, *, /, \, ^, *, /, \ si ^ Semnificatia operatorilor +, -, * este cea de adunare, scadere si înmultire a matricelor cu dimensiuni corespunzatoare Ceilalti operatori au semnificatia din Tabelul 10 Operator Exemplu Semnificatie / a/b a * inv(b) \ a\b inv(a) * b ^ a^b * c=a *b / c = a / b \ c = a \b ^ c = a ^ b Tabelul 10 Operatii cu matrice Prioritatile operatorilor *, /, \ si ^ sunt cele ale operatorilor *, /, \ si ^ In cazul operatorului ^, unul dintre operanzi trebuie sa fie scalar Transpusa unei matrice se calculeaza cu operatorul ' Tabelul 11 arata câteva functii utile pentru lucrul cu vectori si matrice Functia Rezultat length(x) lungimea vectorului x sum(x) suma elementelor vectorului x prod(x) produsul elementelor vectorului x det(A) determinantul matricei A inv(A) inversa matricei A [m n] = size(A) dimensiunile matricei A Tabelul 11 Functii pentru lucrul cu vectori si matrice * în cazul în care matricea este singulara sau aproape singulara, se afisaza un mesaj de eroare * functia size(A) are ca rezultat o matrice cu o linie si doua coloane cu dimensiunile matricei A, m si n Operatorii booleeni ~, &, | si xor se aplica asupra elementelor vectorilor si matricelor In cazul operatorilor &, | si xor matricele trebuie sa aiba aceleasi dimensiuni Exemple Fie matricele : a = [ 0 1 b = si 3 4] Expresia ~b are rezultatul Expresia a & b are rezultatul iar expresia a | b are rezultatul Expresia xor(a, b) are rezultatul Functiile matematice uzuale, sin, cos, etc , pot avea ca argumente vectori si matrice In acest caz functia se aplica asupra fiecarui component al vectorului sau matricei (vectorizarea functiilor) Fie de exemplu matricea A = Avem rezultatele : sin(A) = log(A) = 2 8 2 Generarea vectorilor Se pot genera vectori ale caror elemente cresc cu o valoare data cu operatorul : In prima forma, expresia vi : vf genereaza un vector cu primul element vi si urmatoarele elemente cu valorile vi + 1, vi + 2, , pâna la ultimul element ce are valoare ') afisaza pe ecran prompterul >, citeste o valoare de la tastatura, o atribuie varabilei a si o afisaza pe ecran, ca mai jos a>-2 5 2 12 Functii pentru lucrul cu siruri de caractere Sirurile de caractere în Matlab sunt vectori de caractere Ele se prelucreaza ca si tablourile de numere Caracterele se reprezinta în Matlab în codul ASCII Tabelul 13 prezinta câteva functii utile pentru lucrul cu sirurile de caractere Functia Operatia realizata char(x) converteste un numar într-un caracter num2str (x) converteste numere în siruri de caractere str2num(sir) converteste un sir de caractere într-un numar Tabelul 13 Functii pentru lucrul cu siruri de caractere Expresia x = 'a':1:'f' genereaza vectorul x = cu codurile ASCII ale caracterelor Expresia char(x) converteste numerele în caractere ASCII si genereaza sirul de caractere abcdef Ratia r poate fi si negativa De exemplu, expresia char( ['h':-1:'a']) sau expresia char( 'h':-1:'a') genereaza vectorul hgfedcba Functia poate prelucra si matrice De exemplu, iinstructiunea x = ['e': -1: 'a'; 'a':1:'e'] creaza matricea 101 100 99 98 97 97 98 99 100 101 Expresia char(x) aplicata asupra acestei matrice are rezultatul edcba abcde Functia num2str converteste numere în siruri de caractere Formele functiei sunt num2str(x) si num2str(x, n) unde x este numarul convertit sau matricea convertita în sir de caractere iar n specifica numarul maxim de cifre semnificative convertite Exemplu num2str(2 3134) converteste numarul 2 3134 în sirul '2 3134', în timp ce num2str(2 3134, 3) converteste numarul 2 3134 în sirul '2 31' Exemplu num2str(rand(1, 2), 3) converteste argumentul, matricea rand(1, 2), în sirul de caractere '0 821 0 445' Functia str2num converteste un sir de caractere într-un numar Forma functiei este str2num(sir) unde sir este sirul ce va fi convertit în numar Exemplu Fie instructiunea x = str2num('2 5') + str2num('-3 7') Se convertesc sirurile de caractere în numere si rezultatul este x = -1 2000 2 13 Formatarea rezultatelor Alegerea formatului de afisare a numerelor se face cu instructiunea format cu sintaxa format tip format ('tip') Cele doua forme ale instructiunii sunt echivalente Valorile pentru tip sunt cele din Tabelul 14 Tip Rezultat short Se afisaza 4 cifre dupa virgula, fara exponent long Se afisaza 14 cifre dupa virgula, fara exponent short e Se afisaza 4 cifre dupa virgula, cu exponent long e Se afisaza 14 cifre dupa virgula, cu exponent bank Se afisaza 2 cifre dupa virgula, fara exponent compact Suprima caracterele CR pentru a afisa mai putine linii vide Tabelul 14 Parametrul instructiunii format format short pi/10 ans = 0 3142 format short e pi/10 ans = 3 1416e-001 format bank pi/10 ans = 0 31 format short e format compact pi/10 ans = 3 1416e-001 2 14 Celule Tablouri de celule O celula este un obiect ce poate sa contina obiecte de orice tip, vectori, matrice, alte celule, etc Elementele unei celule se scriu între acolade, { si } De exemplu, o celula ce contine vectorul se defineste astfel : Z = { } Z = Este posibil de a defini tablouri de celule ce pot contine ca elemente celule cu obiecte de tipuri diferite, tablouri, vectori sau matrice, structuri, siruri de caractere, etc Tablourile de celule se definesc asemanator cu tablourile obisnuite, dar cu elementele scrise între acolade, { si } Vectorii de celule sunt utili la scrierea unor mesaje pe mai multe linii Ca exemplu, un vector coloana de celule ce contine ca elemente siruri de caractere este urmatorul A = {'abc' ; 'xyz' ; 'mnp'} Un tablou de celule se creaza cu functia cell() sau prin atribuire De exemplu, instructiunile A = cell(2, 3) B = {[], [], [] ; [], [], []} creaza doua tabele de celule, A si B, cu doua linii si trei coloane cu celule vide Exista doua moduri de a adresa celulele unui tablou : * când utilizam paranteze rotunde, (), adresam celulele tabloului, * când utilizam acolade, {}, adresam valorile celulelor tabloului Fie de exemplu tabloul cu doua celule D = { , [-1+j, 2+3j]} D = Dupa cum se observa, tabloul contine obiecte de tipuri diferite Vrem sa extragem datele din tablou Rezultatul instructiunii C1 = D(2) C1 = este o celula Rezultatul instructiunii C2 = D{2} C2 = -1 0000 + 1 0000i 2 0000 + 3 0000i este chiar data din celula (este un tablou) Vom crea un tablou cu doua celule X, cu aceleasi componente ca D Vom atribui primei celule ca valoare vectorul utilizând paranteze rotunde pentru indice Elementul atribuit trebuie sa fie o celula X(1) = { } X = Vom atribui celei dea doua celule ca valoare vectorul [-1+j, 2+3j] utilizând acolade Elementul atribuit trebuie sa fie un tablou (datele din celula) X{2} = [-1+j, 2+3j] X = Afisarea primului element al vectorului din prima celula a lui X se face cu expresia X{1}(1) X{1}(1) ans = 1 Amintim ca X{1} are ca rezultat vectorul din prima celula, iar X{1}(1) are ca rezultat primul component al acestui vector Fie un vector cu doua celule cu elementele vectorii si [-3 28, 7] Z = { , [-3 28, 7]} Z = Suma celor doi vectori se face ca mai jos Z{1} + Z{2} ans = -2 2800 9 3000 Afisarea elementelor unui tablou de celule se face cu functia celldisp() Ca exemplu, sa afisam vectorul de celule X celldisp(X) X{1} = 1 2 X{2} = -1 0000 + 1 0000i 2 0000 + 3 0000i Pentru a concatena elementele unui vector de celule se utilizeaza operatorul [] Exemplu Concatenarea elementelor vectorului X se face cu expresia [X{:}] ans = 1 0000 2 0000 -1 0000 + 1 0000i 2 0000 + 3 0000i Aplicarea unei functii asupra fiecarei celule a unui tablou se va exemplifica în capitolul referitor la functii 3 Instructiuni Matlab 3 1 Instructiunea for Instructiunea for are forma for ident = vi:r:vf bloc de instructiuni end Aici vi este valoarea initiala a variabilei ident, r este ratia iar vf este valoarea finala Când r > 0, instructiunea for executa blocul de instructiuni repetat cu variabila ident luând valorile vi, vi + r, vi + 2*r, atât timp cât ident = vr Instructiunea for are si urmatoarea forma for ident = expresie bloc de instructiuni end Aici expresie este un vector sau o matrice Când expresie este un vector, ident ia ca valoare câte un element al vectorului si executa blocul de instructiuni Exemplul de mai sus poate fi scris ca j = ; for i = j a(i) = 1 / i ; end Atunci când expresie este o matrice, variabila ident ia pe rând ca valoare câte o coloana a matricei si executa blocul de instructiuni Exemplu Fie matricea Pentru fiecare coloana ai a matricei vrem sa calculam matricea ai * ai' Fie instructiunile a = [ ' [-2 -4 -6]'] for i = a i * i' end Rezultatul este ans = 1 2 3 2 4 6 3 6 9 ans = 4 8 12 8 16 24 12 24 36 Mentionam ca, în prima forma a instructiunii, expresia vi : r : vf este un vector linie Instructiunea for se utilizeaza în acelasi mod în cazul tablourilor de celule Fie vectorul de celule de mai jos D = { , [-1+2j, 2+3j]} D = Afisarea componentele celulelor cu instructiunea for se face astfel for i = 1:2 D{i} end ans = 1 2 ans = -1 0000 + 2 0000i 2 0000 + 3 0000i Reamintim ca D{i} are ca rezultat valoarea celulei 3 2 Instructiunea while Instructiunea while executa repetat un bloc de instructiuni atât timp cât o conditie este adevarata Forma instructiunii este while conditie bloc de instructiuni end Ca exemplu, vom rezolva problema de mai sus i = 1 ; while i = 0 'variabila este pozitiva' else 'variabila este negativa' end 3 5 Instructiunea switch Instructiunea switch executa un bloc de instructiuni în functie de valoarea unei conditii Forma instructiunii switch este switch conditie case expresie bloc case expresie bloc case expresie bloc otherwise bloc end conditia din instructiune este orice expresie aritmetica, relationala sau booleana Pentru început instructiunea switch evalueaza conditia, rezultatul evaluarii este un numar oarecare Când conditia este o expresie booleana sau relationala, rezultatul evaluarii este 0 sau 1 (fals sau adevarat) Instructiunea executa apoi blocul de instructiuni a carui expresie case coincide cu valoarea conditiei Comentariile din programe au forma % comentariu 4 Functii 4 1 Definirea functiilor Functiile se definesc în fisiere cu extensia m Linia de definitie a unei functii este function [y1, y2, , yn ] = numefunctie (x1, x2, xk) Variabilele x1, x2, , xk sunt parametri de intrare Variabilele y1, y2, , yn sunt calculate de functie si sunt parametrii de iesire Parametrii din definitia functiei se numesc formali Functiile pot avea propriile variabile locale Functiile opereaza asupra variabilelor din propriul lor spatiu, adica variabilele locale si parametrii Exemplu de functie ce calculeaza patratul si cubul unei valori function [p, c] = pc(x) c = x * x; p = c * x; O functie este apelata ca termen într-o expresie Parametrii formali din definitia functiei sunt înlocuiti la apelare cu parametrii actuali Parametrii de intrare actuali pot fi expresii ce contin constante, variabile sau functii De exemplu, functia de mai sus este apelata ca [a, b] = pc(2 3) Orice functie are 2 parametri impliciti de intrare: nargin - numarul de parametri de intrare nargout - numarul de parametri de iesire Executia unei functii se termina în mod normal la sfârsitul functiei Pentru a termina executia functiei mai înainte se utilizeaza instructiunea return Functiile sunt salvate în fisiere ce au numele functiei si extensia m Imtr-un fisier m se pot defini mai multe functii, apelate de prima finctie definita în fisier 4 2 Instructiunea global Instructiunea global se utilizeaza pentru a defini variabile globale Daca mai multe functii si spatiul de baza trebuie sa utilizeze aceeasi variabila, acesta variabila trebuie declarata globala, în fiecare functie si în spatiul de baza La atribuirea unei valori acestei variabile în orice functie sau în spatiul de baza, valoarea este accesibila tuturor functiilor si spatiul de baza ce au declarat variabila ca globala Forma instructiunii global este global v1 v2 Exemplu Fie functia urmatoare ce utilizeaza o variabila globala numita comp ce contine o matrice function lm = fcn(x) % declara variabila globala comp global comp lm = x(1) * comp(1, 1) + x(2) * comp(2, 2) ; Programul urmator citeste o matrice în variabila globala comp si apoi apeleaza functia fcn % declara variabila globala comp global comp % citeste datele din fisierul comp txt în matricea comp load comp txt x = ; fcn(x) 4 3 Variabile handle de functii O variabila tip handle de functii contine toate informatiile necesare despre o functie care trebuie apelata mai târziu De obicei, o variabila tip handle de functii este pasata ca argument unei alte functii Instructiunea fh = @numefunctie memoreaza în variabila fh un handle al functiei numefunctie Variabila fh poate fi pasata ca argument altei functii Apelarea functiei numefunctie pasata printr-o variabila handle se face astfel fh (lista de parametri) unde lista de parametri contine parametrii actuali ai functiei numefunctie Exemplu Fie functia fun1(x) care calculeaza valoarea (x + 1) ^ 2 Functia este memorata în fisierul fun1 m function y = fun1(x) y = (x + 1) ^ 2; Sa calculam valoarea fun1(2 3) Functia poate fi apelata în doua moduri: * ca o expresie în membrul drept al unei instructiuni de atribuire, z = fun1(2 3) z = 10 8900 * printr-o variabila handle, fh = @fun1; z = fh(2 3) z = 10 8900 Exemplu Fie functia fun2 ce calculeaza suma si suma patratelor a doua numere Functia este memorata în fisierul fun2 m function [a b] = fun2(x, y) a = x + y; b = x^2 + y^2; Vom apela functia cu valorile x = 2 1, y = -3 3 : * ca o expresie în membrul drept al unei instructiuni de atribuire, [z1 z2] = fun2(2 1, -3 3) z1 = -1 2000 z2 = 15 3000 * printr-o variabila handle, hf = @fun2; [z1 z2] = hf(2 1, -3 3) z1 = -1 2000 z2 = 15 3000 In exemplele anterioare, functiile erau definite în fisiere m Este posibil sa definim în program functii limitate la o singura expresie, fara nume, si sa le atribuim unei variabile handle Ele sunt apelate prin variabila handle, ca mai îninte Asemenea functii fara nume se numesc anonime Definitia unei functii anonime este @(lista parametrilor formali de intrare) expresia functiei Exemplu Sa definim o functie anonima ce calculeaza expresia (x + 1) ^ 2 si sa o apelam cu x = 1 3 fx = @(x)(x+1)^2 ; z = fx(1 3) z = 5 2900 Exemplu Sa definim o functie anonima care calculeaza suma si suma patratelor a doua numere în vectorul [x+y, x ^ 2+y ^ 2] fz = @(x, y) [x+y, x^2+y^2] ; Apelam functia cu valorile x = 2 1, y = -3 3 [z] = fz(2 1, -3 3) z = -1 2000 15 3000 O functie definita într-un fisier m poate avea ca parametru de intrare o variabila handle Variabila handle se poate initializa cu o functie definita într-un fisier m, cu o functie anonima, sau cu o functie standard Matlab In acest fel, o functie poate avea ca parametri alte functii Fie fcn o functie ce are ca parametru o alta functie x = fcn(f, lista de parametri) Aici f este o variabila handle, iar lista de parametri va contine parametrii lui f si alti parametri ai functiei fcn Fie fun o functie ce va fi pasata ca parametru functiei fcn Apelarea functiei fcn se face cu urmatoarea secventa de instructiuni: % creaza variabila handle f f = @fun ; % paseaza variabila handle si % lista de parametri actuali ai functiilor fcn si fun x = fcn(f, lista de parametri actuali); Exemplu Fie functia funb(x) ce va calcula expresia sin(f(x)), unde f(x) este o functie oarecare, ce va fi parametru al functiei funb Functia funb este memorata în fisierul funb m function b = funb(fcn, x) % fcn - variabila handle a unei functii % ce calculeaza o expresie f(x) y = fcn(x); b = sin(y); Vrem sa calculam expresiile sin(x^2) si sin(x^3) cu x = 1 1 Vom defini functia funa ce calculeaza valoarea x^3 ca function a = funa(x) a = x * x * x; si func care calculeaza valoarea x^2 ca function a = func(x) a = x * x; Apelam functia funb cu primul parametru handle al functiei funa si al doilea parametru 1 1 si obtinem fh = @funa z = funb(fh, 1 1) z = 0 9714 Apelam functia funb cu primul parametru handle al functiei func si al doilea parametru 1 1 si obtinem fh = @func z = funb(fh, 1 1) z = 0 9356 Exemplul de mai sus are rolul de a exemplifica utilizarea functiilor Matlab ce au ca parametri alte functii In cazul în care functia ce este parametru este simpla, putem utiliza o functie anonima ca parametru z = funb(@(x)x^3, 1 1) z = 0 9714 Limbajul Matlab are functia standard feval care se utilizeaza pentru evaluarea functiilor standard si a celor utilizator din fisierele m Functia feval utilizeaza o variabila tip handle pentru a evalua functiile Functia feval are forma feval(handle, arg1, arg2, , argn) unde: handle - variabila handle a functiei ce va fi evaluata arg1, arg2, , argn - parametrii functiei ce va fi evaluata Vom rescrie exemplele de mai sus utilizând functia feval Fie functia funb ce calculeaza expresia sin(f(x)) Ea poate avea urmatoarea definitie: function b = funb(fcn, x) % fcn - variabila handle a unei functii % ce calculeaza o expresie f(x) b = sin(feval(fcn, x)); Exemplu Fie functia funa ce calculeaza cubul unei variabile x function a = funa(x) a = x * x * x; Functia funb este apelata pentru a calcula sin(1 13) astfel fh = @funa; b = funb(fh, 1 1); O alta varianta pentru functia funb este urmatoarea function b = funb(fcn, x) % fcn - variabila handle a unei functii % ce calculeaza o expresie f(x) z = feval(fcn, x); b = sin(z); Pentru a calcula sin(x2) este suficient sa definim o functie ce calculeaza patratul unei variabile x si sa o pasam ca parametru functiei funb function a = func(x) a = x * x ; Functia funb este apelata pentru a calcula sin(1 12) astfel: fh = @func; b = funb(fh, 1 1); Functia feval poate evalua si functii anonime De exemplu, pentru a calcula valoarea expresiei x^3 pentru x = 1 5 putem scrie feval(@(x)x^3, 1 5) ans = 3 3750 Matlab are functia predefinita functions ce are ca parametru o variabila tip handle de functii si ca rezultat o structura cu informatiile despre functia continuta în variabila handle De exemplu, pentru a obtine structura cu informatii despre functia funa vom scrie a = functions(@funa) iar rezultatul va fi a = function: 'funa' type: 'simple' file: 'c:\matlab\work\funa' Pentru a afla doar fisierul în care este definita functia vom afisa câmpul file al structurii a în felul urmator a file iar rezultatul va fi ans = c:\matlab\work\funa 4 4 Vectorizarea functiilor Toate functiile standard Matlab pot aplicate asupra vectorilor si matricelor In acest caz functia se aplica asupra fiecarui element al vectorului sau matricei (vectorizarea functiei) Fie de exemplu vectorul a = Instructiunea b = sin(a) are rezultatul 4 5 Aplicarea functiilor asupra tablourilor de celule Functia cellfun() permite sa aplicam o functie asupra fiecarui element al unui tablou de celule Atunci când rezultatul functiei este un scalar si valorile scalare pot fi concatenate într-un tablou, cellfun() are forma A = cellfun(fun, C1, C2, ) unde Ci sunt tablouri de celule, iar fun este o variabila handle de functie ce se plica asupra fiecarie celule din Ci Exemplu Sa calculam lungimea tablourilor din celulele vectorului X cellfun(@length, X) ans = 2 2 Functia cellfun are si forma A = cellfun(fun, C1, C2, 'Param1', val1, ) unde Parami sunt parametri, iar vali sunt valorile acestor parametri Atunci când functia fun nu are ca rezultat un scalar, sau rezultatele nu pot fi concatenate într-un tablou, utilizam parametrul UniformOutput pus la valoarea false Exemplu Sa aplicam functia abs asupra elementelor vectorului de celule X Rezultatul aplicarii functiei abs asupra unei celule este un vector R = cellfun(@abs, X, 'UniformOutput' , false) R = celldisp(R) R{1} = 1 2 R{2} = 1 4142 3 6056 5 Reprezentari grafice 5 1 Reprezentarea curbelor plane Reprezentarea curbelor în plan se face cu functiile din Tabelul 15: Functia Tipul graficului plot scara liniara pe ambele axe loglog scara logaritmica pe ambele axe semilogx scara logaritmica pe axa x semilogy scara logaritmica pe axa y polar coordinate polare stem grafic cu linii terminate cu cercuri stairs grafic în scara bar grafic cu bare Tabelul 15 Functii pentru reprezentarea curbelor plane Functiile au ca parametri perechi de vectori x1, y1, , xn, yn Vectorii x1, x2, , xn contin abscisele curbelor, iar vectorii y1, y2, , yn contin ordonatele curbelor Forma functiei plot este plot(x1, y1, , xn, yn) Proprietatile liniilor cu care se traseaza graficele sunt : culoarea (graficele pot fi trasate cu diverse culori), stilul (continua, etc ) si simbolurile de marcaj (punctele de pe grafic pot fi marcate cu diverse simboluri) Aceste trei proprietati (culoare, simbol de marcaj si stilul liniei) se specifica utilizând simbolurile din Tabelul 16 printr-un sir de caractere, câte unul pentru fiecare proprietate Culoare Simbol Marker Simbol Stil linie Simbol albastru b punct continua - verde g cerc o punctata : rosu r x x întrerupta mov c plus + roz m asterisc * galben y patrat s negru k Tabelul 16 Culori, tipuri de linii si de puncte pentru grafice Pentru a utiliza o anumita culoare, tip de linie sau simbol la curba (xi, yi), perechea de vectori xi, yi este urmata de un sir de caractere din tabelul de mai sus de forma 'cml' ce precizeaza culoarea, markerul si simbolul cu care se traseaza curba De exemplu, functia plot(x, y, 'b') traseaza curba (x, y) cu culoarea albastru Reprezentarea curbelor în mai multe ferestre se face cu instructiunea subplot(m, n, p) Aceasta instructiune împarte ecranul într-o matrice de m*n ferestre si traseaza axele în fereastra p, unde p = 1, , m * n Instructiunea plot urmatoare reprezinta curba în fereastra p Functia xlabel('sir') scrie eticheta axei x (un text sub axa x) Functia ylabel('sir') scrie eticheta axei y (un text dealungul axei y) Functia title('sir') scrie un text ca titlu al figurii Instructiunile precedente pot avea si alti parametrii, care sa modifice tipul, dimensiunea, stilul sau culoarea fontului Un parametru util este FontSize ce da dimensiunea fontului De exemplu, functia xlabel('sir', 'FontSize', 14) scrie eticheta axei x cu dimensiunea fontului 14 Afisarea unei legende pentru cele n curbe dintr-o figura se face cu instructiunile legend('sir1', 'sir2', , 'sirn') sau legend('sir1', 'sir2', , 'sirn', pos) Prima instructiune scrie legenda în coltul dreapta sus A doua instructiune scrie legendele în pozitia specificata de argumentul pos, ce are valorile din Tabelul 17 pos Pozitia 0 cea mai buna pozitie 1 coltul dreapta-sus 2 coltul stânga-sus 3 coltul stânga-jos 4 coltul dreatpa-jos -1 în afara graficului Tabelul 17 Specificarea pozitiei legendei Legenda arata tipul, markerul si culoarea folosite pentru a reprezenta curba, plus textul specificat Mentionam ca, figura, aria de desenare, curba, etichetele axelor, etc , sunt obiecte grafice cu proprietati, create de functiile plot, xlabel, ylabel, etc Toate functiile grafice au ca rezultat un handle al obiectului grafic creat, curba, eticheta axei, legenda, titlul, etc Acest handle poate fi memorat într-o variabila si poate fi utilizat pentru a obtine sau prescrie proprietati ale obiectului grafic De exemplu, pentru a afisa legenda la figuri unde s-a utilizat instructiunea subplot, se utilizeaza variabila handle a instructiunii plot ca mai jos h = plot(t, x) legend(h, 'sir') Functia grid traseaza o grila pe figura Limitele axelor se pot specifica cu instructiunea axis(xmin, xmax, ymin, ymax) Selectarea automata a axelor se face cu instructiunea axis auto Vom exemplifica utilizarea acestor instructiuni pentru reprezentarea grafica a functiilor sin(t) si cos(t) pe intervalul Programul este urmatorul t = 0 :0 001 :2*pi ; plot(t, sin(t), t, cos(t)) xlabel('time') grid title('sin(time), cos(time)') Rezultatul este prezentat în Figura 8 Figura 8 Reprezentarea grafica a functiilor sin si cos pe intervalul Este posibil de a utiliza litere grecesti, alte caractere speciale sau expresii simple în instructiunile ce afisaza siruri: title, xlabel si ylabel Acestea se scriu dupa regulile din Latex Caracterele se definesc în text sub forma \nume caracter Exemple de caractere sunt prezentate în Tabelul 18 Caracter Definitie Caracter Definitie Caracter Definitie \theta \leq \rightarrow ? \Pi ? \neq \uparrow \alpha ' \prime \downarrow ± \pm ? \partial \leftrightarrow \infty \equiv ? \int × \times ÷ \div v \surd + \lfloor + \rfloor \subseteq \supseteq Tabelul 18 Definirea caracterelor speciale Definitia literelor grecesti mari începe cu litera mare, a celor mici cu litera mica, de exemplu \Omega pentru si respectiv \omega pentru Expresiile se definesc dupa urmatoarele reguli: * expresii cu exponenti; expresia ab se defineste ca a^b, * expresii cu indici; expresia xi se defineste ca x i, * daca este cazul, expresiile se grupeaza între acolade, { si } De exemplu, expresia eat se defineste ca e^a^t sau e^{at} Expresia se defineste ca a i iar expresia xi+j se defineste ca x {i+j} Pentru a afisa un text inclinat, (stilul italic), vom defini acel text ca \ittext iar pentru afisare cu litere groase (bold) vom defini textul ca \bftext Pentru exemplificare vom reprezenta grafic functia sin(?) pe intervalul Vom scrie urmatoarele texte: * pe axa y, textul sin(?), * pe axa x, textul 0 = ? = 2 p, * titlul graficului va fi Functia sin(?) Programul este urmatorul t = 0: 0 001: 2*pi; x = sin(t); plot(t, x) xlabel('0 \leq \theta \leq 2\pi') ylabel('sin(\theta )') title('Functia sin(\theta)') grid Rezultatul programului este cel din Figura 9 Figura 9 Reprezentarea grafica a functiei sin(?) pe intervalul Pentru a utiliza interpreterul Latex generam un sir de caractere cu expresia ce dorim sa fie afisata si apelam functia text Pentru exemplificare vom reprezenta functia sin(t)/t Legenda axei x, t[s], va fi scrisa cu litere groase (bold) xlabel('\bft[s]') Vom afisa pe figura formula sin(t) / t în centrul figurii Pentru aceasta vom crea un sir de caractere cu formula scrisa în Latex texstr = '$\frac{sin(t)}{t}$'; Sirurile de caractere se afisaza pe grafic într-o anumita pozitie cu functia text cu forma text('NumeProprietate', ValoareProprietate, ) Proprietatile des utilizate ale obiectului text sunt : * 'String' - sirul ce va fi afisat, * 'FontSize' - da dimensiunea fontului în unitati de font, * 'FontUnits' - unitatile de font, 'points', 'normalized', 'pixels', 'centimeters', * 'HorizontalAlignment' - aliniamentul orizontal, 'left', 'center', right, * 'VerticalAlignment' - aliniamentul vertical, 'top', 'middle', 'bottom', * 'Interpreter' - interpreterul utilizat, 'latex', 'tex', 'none', * 'Units' - unitatile de masura, 'points', 'normalized', 'pixels', 'centimeters' * 'Position' - pozitia sirului, [x, y] Optiunea 'normalized' a proprietatii 'Units' transforma coltul din stânga jos al dreptunghiului figurii în (0, 0) si cel din dreapta sus în (1, 1) Vom apela functia text astfel : text('string', texstr, 'interpreter', 'latex', 'fontsize',40, 'units','norm', 'pos',[ 5 5]); Programul este urmatorul t = 0:0 01:4*pi; x = sin(t) / t; plot(t, x) grid xlabel('\bft[s]') texstr = '$\frac{sin(t)}{t}$'; text('string',texstr, 'interpreter', 'latex', 'fontsize',40, 'units','norm', 'pos',[ 5 5]); Rezultatul programului este cel din Figura 10 Figura 10 Reprezentarea grafica a functiei sin(t) / t Implicit, fiecare instructiune plot creaza o noua figura Desi o instructiune plot poate reprezenta grafic mai multe curbe, uneori este necesar sa reprezentam mai multe curbe pe aceeasi figura cu mai multe instructiuni plot Instructiunea hold on mentine figura curenta si proprietatile axelor, astfel încât urmatoarele instructiuni plot reprezinta curbe pe aceeasi figura Instructiunea hold off reface modul implicit al instructiunii plot, de a sterge graficele precedente si de a reface proprietatile initiale ale axelor înainte de a desena un noi grafice Exemplu Fie programul: t = 0: 0 01: pi; a = sin(2 * t); b = cos(4 * t); plot(t, a); hold on plot(t, b) El traseaza doua curbe cu doua instructiuni plot pe aceeasi figura Axele pot fi facute vizibile sau invizibile cu optiunile on si off: axis on axis off In unele cazuri avem nevoie sa reprezentam simultan mai multe figuri Dupa cum am spus mai înainte, acest lucru se face cu instructiunea subplot (m, n, p) care împarte figura într-o matrice de mxn axe si selecteaza axa p pentru reprezentarea grafica ce se va face cu instructiunea plot urmatoare Ca exemplu vom reprezenta curbele sin(t), cos(t), sin2(t) si cos(t) în functie de sin(t) într-o matrice de patru figuri Porgramul este urmatorul t = 0: pi/20: 2*pi; x = sin(t); y = cos(t); z = x * x; subplot(2, 2, 1) plot(t, x) grid xlabel('t[s]') ylabel('sin(t)') subplot(2, 2, 2) plot(t, y) grid xlabel('t[s]') ylabel('cos(t)') subplot(2, 2, 3) plot(t, z) grid xlabel('t[s]') ylabel('sin^2(t)') subplot(2, 2, 4) plot(x, y) grid xlabel('sin(t)') ylabel('cos(t)') Rezultatul rularii programului este cel din Figura 11 Figura 11 Reprezentare grafica într-o matrice de figuri 5 2 Editarea unei figuri Figurile si componentele lor, aria de desenare, curbele, etichetele axelor, titlul, etc , sunt obiecte cu proprietati ce pot fi modificate Modificarea proprietatilor obiectelor se poate face în doua feluri : * interactiv, cu editorul de proprietati, * prin program, cu instructiuni ce modifica proprietatile Vom prezenta în continuare cele doua moduri de prelucrare a figurilor Bara de meniu a unei figuri este cea din în Figura 12 Salvarea unei figuri se face cu optiunile meniului File Editarea figurii se face cu optiunile meniului Edit Inserarea de text, sageti sau linii, titlu, etc se poate face cu optiunile meniului Insert sau cu butoanele barei de instrumente a figurii Optiunile acestor meniuri sunt aratate în Tabelul 20 Figura 12 Bara de meniu si bara de instrumente ale unei figuri a) meniul File b) meniul Edit c) meniul Insert Tabelul 19 Meniurile File, Edit si Insert ale unei figuri Salvarea unei figuri se face cu optiunile Save sau Save As ale meniului File La editarea unei figuri se pot adauga text, linii sau sageti Pentru aceasta se utilizeaza butoanele barei de instrumente din Figura 12, vezi si Figura 13, sau optiunile Arrow, Line si Text ale meniului Insert din Tabelul 19 Figura 13 Butoane pentru editarea unei figure Semnificatia acestor butoane si a optiunilor meniului Insert este cea din Tabelul 20 Buton Meniul Insert Semnificatie Selecteaza figura Text Adauga text Arrow Adauga o sageata Line Adauga o linie Tabelul 20 Inserarea de texte, linii sau sageti Procedura de adaugare a unui text, a unei sageti sau a unei linii cu butoanele barei de instrumente este urmatoarea : 1 se selecteaza figura cu un clic pe butonul sau prin afisarea meniul Insert, 2 clic pe unul din butoanele , sau , sau pe optiunile Text, Arrow, sau Line ale meniului Insert, 3 se introduce textul, sageata sau linia dorita 5 3 Modificarea interactiva a proprietatilor obiectelor grafice O reprezentare grafica este formata din obiecte grafice ce formeaza un arbore La primul nivel se afla un obiect numit root care reprezinta ecranul Pe urmatorul nivel se gasesc obiectele tip figure (figurile) Obiectele tip figure au ca si componente obiecte tip axes care reprezinta axele (zona de desenare) Obiectele tip axes au ca si componente obiecte tip line ce reprezinta curbele Acest arbore este afisat de Editorul de proprietati ca în Figura 14 Figura 14 Arborele unei figuri afisat de Editorul de proprietati Obiectele interesante sunt axes si line ce permit modificarea proprietatilor axelor si curbelor 5 3 1 Modificarea proprietatilor axelor Obiectul axes reprezinta zona de desenare : axele, etichetele axelor si titlul figurii, etc Proprietatile obiectului axes sunt : tipul (numele fontului), dimensiunea si stilul fontului utilizat în afisarea textelor de pe axe si din zona de desenare, unitatea de masura (pixeli, inci, centimetri, etc ), limitele axelor, pozitia axelor, culoarea zonei de desenare, etc Modificarea acestor proprietati se face cu editorul de proprietati Editorul de proprietati se lanseaza în mai multe moduri : * clic pe optiunea Axes Properties a meniului Edit, * clic pe optiunea Figure Properties a meniului Edit si în caseta Edit Properties for se alege optiunea axes, vezi Figura 14, * se afisaza meniul contextual din Figura 15 dupa procedura de mai jos si clic pe optiunea Properties 1 se selecteaza figura cu un clic pe butonul , 2 se executa un clic cu butonul stâng pe axe sau pe figura pentru selectare, 3 se executa un clic cu butonul drept pe figura, curba, legendele axelor, titlul figurii care afisaza un meniu contextual In Figura 15 este afisat meniul contextual afisat la un clic pe figura 4 se executa un clic pe optiunea Properties a acestor meniuri, ce afisaza caseta editorului de proprietati Figura 15 Meniu contextual pentru editarea unei figuri Caseta editorului de proprietati pentru modificarea axelor, Property Editor - Axes, are câteva sectiuni importante, sectiunea Style din Figura 16 si sectiunile X, Y, Z, vezi Figura 17 pentru sectiunea X Modificarea stilului (proprietatilor) axelor, font (dimensiune, tip), etc , se face în sectiunea Style a casetei de dialog Property Editor - Axes din Figura 16 In aceasta sectiune se pot modifica : * titlul figurii în caseta de text Title, * stilul textului titlului se alege în sectiunea text a editorului afisata cu un clic pe butonul Properties, * culoarea fondului se alege din lista Background, * grosimea liniilor axelor se alege din lista Axes line width, * fontul, tip si dimensiune pentru axe se aleg din casetele Font name, Font size, Font units, Font angle, etc Figura 16 Caseta Property Editor, sectiunea Style Modificarea limitelor axelor se face în sectiunile X, Y, Z ale casetei de dialog corespunzatoare Caseta axei X este aratata în Figura 17 Atunci când optiunea Limits este validata, limitele sunt alese automat astfel încât este reprezentata întreaga figura Pentru a modifica limitele pentru a afisa doar o parte din figura se invalideaza optiunea Limits si se aleg limitele dorite Se mai pot modifica : * eticheta axei în caseta de text Label, * stilul textului etichetei axei, tipul fontului, dimensiune, culoare, etc , se alege în sectiunea Text a editorului afisata cu un clic pe butonul Properties, * culoarea axei se alege din lista Color, * scala, lineara sau logaritmica, se alege validând unul din butoanele din zona Scale, * pozitia axei se alege din lista Location Figura 17 Caseta Property Editor, modificarea limitelor axelor 5 3 2 Modificarea proprietatilor curbelor Modificarea proprietatilor liniei cu care se traseaza curba se face cu caseta de text Property Editor - Line din Figura 18 In sectiunea Style a acestei casete se pot alege: 1 în sectiunea Line Properties: * stilul liniei din lista Line style, * grosimea liniei din lista Line width, * culoarea liniei din lista Line color 2 în sectiunea Marker Properties: * simbolul cu care se marcheaza punctele curbei din lista Style, * dimensiunea simbolului din lista Size, * culoarea simbolului din lista Edge color Figura 18 Caseta Property Editor pentru modificarea proprietatilor curbei Mentionam ca o serie de proprietati ale figurii, axelor, liniei se pot alege din optiunile meniurilor contextuale afisate la un clic cu butonul drept pe figura, scala axelor, etichetele axelor Aceste meniuri sunt aratate în Tabelul 21 a) clic pe eticheta unei axe sau titlul figurii b) clic pe figura sau axele figurii c) clic pe curba, text, linie, sageata Tabelul 21 Meniuri contextuale afisate la un clic pe figura In toate cazurile, un clic pe optiunea Properties afisaza caseta editorului de proprietati In meniul din Tabelul 21a se pot modifica dimensiunea, stilul si culoarea fontului In meniul din Tabelul 21c se modifica grosimea, stilul si culoarea curbei, a unei sageti, linii sau text Optiunile Cut, Copy, Paste efectueaza operatiile respective cu Clipboardul Optiunea Clear sterge obiectul respective In Figura 19 este aratata o figura editata S-au modificat stilul si dimensiunea fontului axelor, eticheta axei X si titlul figurii Figura 19 Exemplu de figura editata Un alt exemplu de figura editata este cea din Figura 20 Aici s-au inserat o sageata si un text si s-au modificat stilul acestora si al curbei Inserarea obiectelor s-a facut utilizând bara de instrumente a figurii Modificarea proprietatilor se poate face cu editorul de proprietati s-au utilizând optiunile meniurilor contextuale afisate cu un clic cu butonul drept al mousului pe obiect Figura 20 Exemplu de figura editata In cazul utilizarii Editorului de proprietati se afisaza arborele figurii, vezi Figura 21, apoi se afisaza casete de dialog corespunzatoare obiectelor de unde se aleg proprietatile dorite Figura 21 Ierarhia obiectelor grafice ale unei figuri 5 4 Modificarea programata a proprietatilor obiectelor grafice O figura este compusa din mai multe obiecte grafice Obiectele grafice formeaza un arbore care are ca radacina un obiect numit Root care corespunde ecranului Figurile, (obiecte tip Figure), sunt urmatoarele obiecte grafice De regula, ele sunt create automat de functiile pentru reprezentari grafice, plot, etc Toate figurile sunt copiii obiectului Root Orice alte obiecte grafice sunt copii ai unei figuri Obiecte interesante ale unei figuri sunt Axes, Line si Text Obiectele tip Axes definesc axele, o regiune în figura unde sunt reprezentate curbele Obiectele Axes sunt parinti ai obiectelor de tip Line si Text Obiectele tip Line sunt create de functiile pentru reprezentari grafice, plot, etc si reprezinta curbele Obiectele tip Text sunt create de functiile title, xlabel si ylabel reprezinta siruri de caractere Fiecare obiect grafic este identificat de o variabila handle care este rezultatul functiei ce îl creaza Utilizând aceasta variabila handle se pot modifica proprietatilor obiectului grafic Pentru modificarea proprietatilor unui obiect se utilizeaza functia set cu forma set(handle-obiect, 'NumeProprietate', 'ValoareProprietate') Pentru a obtine valoarea unei proprietati a unui obiect se utilizeaza functia get cu forma get(handle-obiect, 'NumeProprietate') In aceste instructiuni handle-obiect este variabila handle a obiectului Pentru a obtine valoarea variabilei handle-obiect pentru diverse obiecte se utilizeaza functiile din Tabelul 22 Obiect grafic Functie Figure gcf Axes gca Obiectul curent selectat gco Tabelul 22 Functii ce dau variabila handle pentru obiecte grafice Instructiunea get(handle-obiect) afisaza toate proprietatile obiectului identificat de variabila handle-obiect Ierarhia obiectelor grafice este cea din Figura 22 Dupa cum se vede din Figura 22 obiectele de tip Line si Text sunt copiii obiectului Axes Figura 22 Ierarhia obiectelor grafice Aceasta ierarhie este afisata si de Editorul de proprietati, vezi Figura 14 In continuare vom prezenta proprietatile importante ale obiectelor grafice Toate obiectele grafice au proprietatile comune din Tabelul 23 Proprietate Continut Parent Variabila handle-obiect a parintelui Children Vector cu variabile handle-obiect ale copiilor Type Tipul obiectului Tabelul 23 Proprietati comune obiectelor grafice Stergerea unui obiect grafic se face cu functia delete(handle-obiect) 5 4 1 Obiectul Root Obiectul Root reprezinta ecranul Proprietatile interesante ale acestui obiect sunt cele din Tabelul 24 Proprietate Valori ScreenSize [left, bottom, width, height] Units {pixels, inches, centimeters, points} Tabelul 24 Proprietatile obiectului Root Obiectul Root are handle 0 Toate celelalte obiectele grafice au proprietatile comune referitoare la font din Tabelul 25 Proprietate Valori FontAngle {normal, italic, oblic} FontName {Courier, FixedWidth} FontSize Dimensiunea fontului în FontUnits FontUnits {points, normalized, pixels, inches} FontWeight {normal, bold, light} Tabelul 25 Proprietati commune obiectelor grafice relative la font 5 4 2 Obiectul Figure Proprietatile obiectului Figure sunt cele din Tabelul 26 Proprietate Valoare Position [left, bottom, width, height] Units {pixels, inches, centimeters, points} DefaultAxesColor {b, g, r, c, m, y, k} CurrentAxes Sistemul curent de axe Tabelul 26 Proprietatile obiectului Figure Proprietatea Position este un vector cu patru componente ce dau pozitia fata de coltul stânga jos al ecranului si dimensiunile figurii Unitatea de masura este data de proprietatea Units Sistemul curent de axe este variabila handle a obiectului Axes creat cu functia axes De regula, figurile sunt create de functiile pentru reprezentari grafice, plot, etc Obiectul figure poate fi creat cu functia hf = figure ('proprietate', 'valoare', ) unde hf este variabila handle a obiectului La crearea unei figuri ne intereseaza proprietatea ScreenSize a obiectului Root Ea se obtine ca scrsz = get(0, 'ScreenSize) rezultatul este un vector cu patru componente, [left, bottom, width, height] La crearea unei figuri, aceste componente pot fi utilizate pentru a defini pozitia figurii pe ecran 5 4 3 Obiectul Axes Obiectul Axes corespunde zonei de desenare în care sunt reprezentate axele, curbele si textele afisate pe axe sau în zona de desenare Obiectul axes se creaza cu instructiunea haxes = axes('proprietate', 'valoare', ) unde haxes este variabila handle a obiectului Orice figura are un sistem curent de axe (aria curenta de desenare) Acest sistem de axe este proprietatea CurrentAxes a figurii si este variabila handle a unui sistem de axe creat cu functia axes De exemplu, instructiunile ha = axes( ) set(gcf, 'currentaxes', ha) creaza un sistem de axe si îl face sistemul curent de axe al figurii Obiectul Axes are proprietatile din Tabelul 27 Ele reprezinta pozitia zonei de desenare în figura, pozitia axelor, etichetele axelor si titlul, culoarea axelor si a textului Proprietate Valoare Position [left, bottom, width, height] OuterPosition [left, bottom, width, height] Units {pixels, inches, centimeters, points} XAxisLocation {top, bottom} YAxisLocation {left, right} Color {none, [r, g, b]} XColor {b, g, r, c, m, y, k} YColor {b, g, r, c, m, y, k} XLabel Variabila handle-obiect a axei X YLabel Variabila handle-obiect a axei Y Title Variabila handle-obiect a titlul figurii Tabelul 27 Proprietatile obiectului Axes Proprietatea Position este un vector cu patru componente ce specifica dreptunghiul ce încadreaza zona de desenare left si bottom dau pozitia coltului stânga jos al dreptunghiului fata de coltul stânga jos al figurii width si height dau dimensiunile dreptunghiului, vezi Figura 23 Unitatea de masura este data de proprietatea Units Proprietatea OuterPosition este tot un vector cu patru componente ce specifica dreptunghiul ce incadreaza zona de desenare, axele si titlul left si bottom dau pozitia coltului stânga jos al dreptunghiului fata de coltul stânga jos al figurii width si height dau dimensiunile dreptunghiului Unitatea de masura este data de proprietatea Units Pentru a modifica, de exemplu, dimensiunea fontului axelor la valoarea 18 se utilizeaza functia set(gca, 'FontSize', 18) Figura 23 Pozitia axelor în cadrul unei figuri Celelalte proprietati ale obiectului Axes sunt evidente Codurile culorilor sunt cele din Tabelul 16 Exemple de utilizare a acestor proprietati vor fi prezentate ulterior 5 4 4 Obiectul Line Obiectul Line reprezinta o curba Proprietatile acestui obiect sunt cele din Tabelul 28 Proprietate Valoare XData Vector cu abscise YData Vector cu ordonate Color {b, g, r, c, m, y, k} Tabelul 28 Proprietatile obiectului Line Un obiect Line este creat cu instructiunea line de forma hline = line ('proprietate', 'valoare', ) unde hline este variabile handle a obiectului 5 4 5 Obiectul Text Obiectul Text reprezinta un sir de caractere creat si afisat de instructiunile xlabel, ylabel, title Obiectul Text are proprietatile din Tabelul 29 Proprietate Valoare String 'text' HorizontalAlignment {left, center, right} VerticalAlignment {top, middle, bottom} Interpreter {latex, tex, none} Units {points, normalized, pixels, centimeters} Position [x, y] Tabelul 29 Proprietatile obiectului Text Optiunea 'normalized' a proprietatii 'Units' transforma coltul din stânga jos al dreptunghiului figurii în (0, 0) si cel din dreapta sus în (1, 1) Un obiect text este creat cu functia text are forma ht = text(x, y, 'sir de caractere') sau forma ht = text(x, y, 'sir de caractere', 'NumeProprietate', 'ValoareProprietate') Functia text are ca rezultat variabila handle al obiectului creat, ht Functia scrie textul dintre apostrofuri în punctul de coordonate (x, y) al figurii x si y pot fi vectori, caz în care textul este scris în toate perechile de coordonate date de x si y Daca textul este un tablou cu acelasi numar de linii ca dimensiunea vectorilor x si y, atunci functia text afisaza în fiecare punct linia corespunzatoare din tablou Parametri NumeProprietate si ValoareProprietate dau caracteristicile textului: tipul si dimensiunea fontului, etc Exemplu Vom reprezenta grafic functia sin(t) creind direct figura, aria de desenare, si curba Vom crea o figura cu latimea jumatatea din latimea ecranului si înaltimea jumatatea din înaltimea ecranului Pozitia punctului din stânga jos va fi (1, latimea ecranului / 2) Dimensiunile ecranului, în pixeli, se obtin cu functia scrsz = get(0, 'screensize') Figura se creaza cu instructiunea hf = figure('position', ) Vom crea o arie de desenare cu instructiunea ha = axes('position', ) pentru a avea loc pentru etichetele axelor Vom reprezenta grafic o functie cu instructiunea hl = line(t, st) Programul este cel de mai jos t = 0: 0 01: 2*pi; st = sin(t); scrsz = get(0, 'screensize') scrsz = 1 1 1440 900 hf = figure('position', ) hf = 1 ha = axes('position', ) ha = 158 0076 hl = line(t, st) hl = 159 0076 thx = text('string', 't') thx = 160 0076 set(ha, 'xlabel', thx) thy = text('string', 'sin(t)') thy = 161 0077 set(ha, 'ylabel', thy) grid Rezultatul fiecarei functii a programului este cel din figurile de mai jos Mentionam ca, odata curba creata, etichetele axelor se puteau crea si cu functiile xlabel si ylabel Figura 24 Figura creata cu instructiunea figure Figura 25 Figura si aria de desenare Figura 26 Figura, aria de desenare si curba Figura 27 Figura, aria de desenare, curba si etichetele axelor 5 4 6 Modificarea proprietatilor obiectelor grafice In acest paragraf vom exemplifica modificarea proprietatilor obiectelor grafice ce constituie o figura Putem de exemplu afla eticheta axei X astfel : % fa este variabila handle-obiect a obiectului axes fa = get(gca) % afisaza proprietatile obiectului axes get(fa) % fxlabel este variabila handle-obiect a obiectului axa X fxlabel = get(fa, 'xlabel') % obtine valoarea proprietatii 'string' a obiectului axa X % care este eticheta axei X get(fxlabel, 'string') Vom exemplifica utilizarea obiectelor grafice cu un program care reprezinta grafic functia sin(?) pe intervalul Vom scrie urmatoarele texte: * pe axa y textul sin(T), * pe axa x textul 0 = ? = 2 p, * textul titlului va fi Plot of sin(T), * în punctul (p, 0) vom afisa mesajul 0 Functia 1(t) este functia treapta unitara Ea este predefinita ca heaviside(t) syms t f F syms a positive f = exp(-a*t)*heaviside(t) f = exp(-a*t)*heaviside(t) F = fourier(f) F = 1/(a+i*w) Pentru verificare vom calcula transformata Fourier inversa a functiei F ifourier(F) ans = exp(-a*x)*heaviside(x) 19 2 Transformata Laplace Transformata Laplace directa se defineste cu formula unde f(t) este o functie simbolica scalara de argument t, iar F(s) este transformata Laplace de argument s Transformata Laplace directa simbolica a functiei f se calculeaza cu functia: F = laplace(f) F = laplace (f, v) F = laplace (f, u, v) In prima forma, variabila functiei f este x, iar a functiei F este w (variabilele implicite) In forma a doua, variabila functiei f este t, (variabila implicita), iar a functiei F este v In forma a treia, variabila functiei f este u, iar a functiei F este v Transformata Laplace inversa se defineste cu formula unde F(s) este o functie simbolica scalara de argument s, iar f(x) este transformata Laplace inversa de argument t Transformata Laplace inversa simbolica a functiei F se calculeaza cu functia: f = ilaplace(F) f = ilaplace (F, u) f = ilaplace (F, v, u) In prima forma, variabila functiei F este s, iar a functiei f este t (variabilele implicite) In forma a doua, variabila functiei F este s, (variabila implicita), iar a functiei f este u In forma a treia, variabila functiei F este v, iar a functiei f este u Exemplu Sa calculam transformata Laplace a functiei syms f F t s syms a positive f = exp(-a *t) f = exp(-a*t) F = laplace(f) F = 1/(s+a) Exemplu Sa calculam transformata Laplace inversa a functiei syms f F s t syms a positive F = 1/(s+a) f = ilaplace(F) f = exp(-t*a) 19 3 Transformata Z Transformata Z directa se defineste cu formula unde f este o functie simbolica scalara de argument n, iar F(z) este transformata Z de argument z Transformata Z directa se calculeaza cu functia: F = ztrans(f) F = ztrans(f, w) F = ztrans(f, k, w) In prima forma variabila functiei f este n, iar a functiei F este z, In a doua forma variabila functie f este n (variabila implicita), iar variabila functiei F este w In a treia forma variabila functiei f este k, iar a functiei F este w Transformata Z inversa se calculeaza cu formula unde F(z) este o functie simbolica scalara de argument z, iar f(n) este transformata Z inversa de argument n Transformata Z inversa se calculeaza cu functia: f = iztrans(F) f = iztrans(F, k) f = iztrans(F, w, k) In prima forma, variabila functiei F este z, iar a functiei f este n (variabilele implicite) In forma a doua, variabila functiei F este z, (variabila implicita), iar a functiei f este k In forma a treia, variabila functiei F este w, iar a functiei f este k In exemplele urmatoare vom nota cu T perioada de esantionare Momentele de esantionare vor fi nT Exemplu Sa calculam transformata Z a functiei Functia de timp ce este esantionata este f(t) = e-at syms a n T f = exp(-a*T*n) f = exp(-a*T*n) F = ztrans(f) F = z/exp(-a*T)/(z/exp(-a*T)-1) Vom calcula acum transformata Z inversa a functiei F iztrans(F) ans = exp(-a*T)^n Exemplu Sa calculam transformata Z inversa a functiei syms a z; F = z / (z - a); f = iztrans(F) f = a^n Functia esantionata este f(n) = an Cu notatia obtinem In consecinta, functia de timp ce este esantionata este , unde Exemplu Sa calculam transformata Z inversa a functiei syms a z; F = z / (z + a); f = iztrans(F) f = (-a)^n Exemplu Sa se calculeze transformata Z inversa a functiei syms z F F = z/(z-1)^2 F = z/(z-1)^2 f = iztrans(F) n Functia de timp ce este esantionata este f(t) = t/T Vom calcula acum transformata Z a functiei f ztrans(f) ans = z/(z-1)^2 33 